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Abstract 

Diamond relay channel model, as a basic transmission model, has recently been attracting consider- 
able attention in wireless Ad Hoc networks. Node cooperation and opportunistic scheduling scheme are 
two important techniques to improve the performance in wireless scenarios. In the paper we consider 
C*~) . such a problem how to efficiently combine opportunistic scheduling and cooperative modes in the 

r» ■ 

Rayleigh fading scenarios. To do so, we first compare the throughput of SRP (Spatial Reused Pattern) 
C ' and AFP (Amplify Forwarding Pattern) in the half-duplex case with the assumption that channel side 

o ; 

Q\ t information is known to all and then come up with a new scheduling scheme. It will that that only 

switching between SRP and AFP simply does little help to obtain an expected improvement because 
SRP is always superior to AFP on average due to its efficient spatial reuse. To improve the throughput 
' further, we put forward a new processing strategy in which buffers are employed at both relays in 

SRP mode. By efficiently utilizing the links with relatively higher gains, the throughput can be greatly 
improved at a cost of queuing delay. Furthermore, we shall quantitatively evaluate the queuing delay and 
the tradeoff between the throughput and the additional queuing delay. Finally, to realize our developed 
strategy and make sure it always run at stable status, we present two criteria and an algorithm on the 
selection and adjustment of the switching thresholds . 

Index Terms 

diamond relay networks, cooperative pattern, block Rayleigh fading, finite state channel, G/G/l 
queueing system, tradeoff. 
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I. Introduction 

In recent years, motivated by higher quality-demanded applications in the wireless Ad Hoc 
networks, cooperation among nodes is considered to be more useful over various relay network 
models. Early in 1970s, a classical three-node relay channel comprised of one source, one 
destination and only one relay was first introduced by van de Meulen in [3]. Then, in reference 
[4] and [5], the cut-set bound and the achievable rates with power allocation were studied for the 
half-duplex case in this three-node network. Especially, Reference [5] presented some lower and 
upper bounds on outage capacity of three-node network. The analysis of diversity-multiplexing 
tradeoff for the three-node network in the half-duplex case was given in [6]. Besides, systems 
using multiple relays were also studied. In [7], both the achievable rates and the upper bound 
of capacity were studied for the diamond relay network in the half-duplex case. The two-relay 
model (i.e. the Diamond Relay Networks) operating in the full-duplex case were studied in 
[8]. In [9], they analyzed a model in which cooperative communication proceeds in a parallel 
relay network where the exogenous arrival of packets and the FIFO (First In First Out) queueing 
system are introduced. Multiple relays using orthogonal channels were analyzed in [10]. Besides, 
relay networks with and without delay were discussed in [12] and [13], respectively. Reference 
[14]-[19] presented analysis of capacity and delay tradeoff in the networks comprised of many 
random prelocated nodes. These works showed that a cooperative gains can be obtained in 
distributed wireless networks if nodes can help each other to relay information. This motivated 
us to study the wireless relay networks further, especially for some classical topologies and 
cooperative modes. 

Referring to the wireless relay networks, the three-node network has been a very hot topic in 
the research area of the cooperative network since 1970s due to its classical and representative 
topology. In the late ninety's, B.Schein and R.Gallager proposed another kind of relay network, 
the Gaussian parallel relay network [8], in which the diamond relay network was first introduced 
implicitly. The reasons that diamond relay networks were considered include two facets: Firstly, 
this kind of topology is relatively easy traceable in theory and has more freedom than three-node 
network model. Secondly, this model can be used in some wireless scenarios, where a sender is 
convenient to select a few neighbors. Fig. 1 shows one application scenario where the black and 
square parts represent the obstacles blocking radio signals seriously, e.g., skyscrapers in business 
district. In such a scenarios, if node A wants to transmit to node D at a higher rate, the diamond 
relay network comprised of node A,B,C and D can do a great help. 

Based on different topologies of the wireless relay networks [3]-[8],[10], several cooperative 
modes were introduced consequently, such as SRP (Spatial Reused Pattern) and AFP (Amplify 
Forwarding Pattern) in [7]. In the wireless time-varying and fading scenarios, different coopera- 
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tive modes may have big difference in term of reliable throughput where opportunistic scheduling 
may play a key role in improving the throughput. Motivated by this, we firstly review the two 
different cooperative modes, namely SRP and AFP, and present their corresponding capacities in 
the diamond relay network, shown in Fig. 2, in the half-duplex case. In the discussion, relay nodes 
adopt two relay schemes, namely, the amplify-and-forward (AF - relay node simply amplifies the 
signals received from source and forwards to destination) and decode-and-forward (DF - relay 
node decodes the information received from source, re-encodes and forwards to destination). In 
the sequel, AF relay scheme is referred to AFP (Amplify Forwarding Pattern) mode while DF 
relay scheme is referred to SRP(Spatial Reused Pattern) mode. Furthermore, we compare the 
throughput of the two cooperative modes and find that there exists a big difference between the 
performances of the SRP and AFP mode. In some cases, SRP has a larger throughput than AFP. 
Otherwise, it has smaller throughput than AFP. Thus, our first contribution is that we analyze 
the efficiency of opportunistic scheduling and put forward a hybrid relay scheme with switching 
between SRP and AFP so that it can be adapted to the channel variation. 

Although the opportunistic scheduling is usually considered to be effective, it will be shown 
that combining SRP and AFP in a simple way can not obtain an expected improvement. This is 
because SRP is always superior to AFP on average due to the efficient spatial reuse. In previous 
works [3]-[8], using Max-Flow Min-Cut theorem, it can be seen that the whole performance 
is always reduced greatly by some bottleneck links usually caused by the fading and time- 
vary characteristics of the wireless channels. Therefore, using buffers at nodes maybe help to 
efficiently employ some channels in rather good conditions. In this way, it will improve the whole 
performance of the throughput in the wireless networks at the cost of some additional queuing 
delay. Motivated by this, we shall put forward a new processing strategy in which buffers are 
employed at both relays in SRP mode. Its basic idea is that in diamond relay channels, sometimes 
there is only one or two links in very good conditions while the others in bad ones. In this case, 
buffers can help relays efficiently utilizing the links in good conditions so that the throughput 
is greatly improved at some cost of queuing delay. One important contribution of this work is 
that we quantitatively evaluate the queuing delay and discuss its tradeoff with the throughput 
where two criteria are considered to characterize the thresholds for separating pretty good or 
rather bad condition of channels and make sure the network run in stable status. In addition, an 
adjust algorithm is also given. 

For simplicity, finite state channel model is used in which the received SNR is partitioned into 
N levels. The status of the four links are assumed to be independent and each link operates in a 
memoryless mode. G/G/l queueing system is used to model the relay node with buffers and the 
associated two links connecting to it. Besides, a continuous traffic model is considered in which 
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the source always has information to send. Then one can use Marshall's theory to solve it. Note 
that here we consider the source always has information to send, which will help us to get more 
insights on the maximum achievable throughput of this diamond network. Thus, other traffic 
models, such as Bernoulli or Poisson traffic model, will not be discussed here. Another point 
should be mentioned is that due to the IC design improvement, buffers with enough capacity 
are becoming much cheap and with low cost, so the delay is a more rigorous factor influencing 
the performance of the networks. Therefore we emphasize the average queuing time without 
limitation of the buffer sizes, referred to delay in this paper. 

The rest of this paper is organized as follows: In Section HH we introduce the system model. 
In Section Unl we firstly compare the performance of SRP and AFP scheme,then we give an 
opportunistic scheduling scheme, namely an hybrid scheme of combining SRP and AFP modes. 
In Section [TV] we shall propose a new processing strategy in which buffers are employed at the 
relay nodes in SRP mode. We present two criteria to characterize the thresholds for realizing 
tradeoff between the throughput and queuing delay and make sure the network run in stable 
status. Some simulation results are given in Section |V] Finally, We present the conclusion in 
section |VT1 

II. System Model 

A. Diamond Relay Channel and Two Time Sharing Patterns: SRP and AFP 

The discussion of diamond relay channel was first considered in Schein and Gallager's work 
in a full-duplex case [8]. Recently, Feng, etc, studied this kind of networks in the half-duplex 
case in [7], shown in Fig. 2. It is comprised of four nodes including a source N s , two relays 
Ni and N 2 and a destination N d . It is assumed that all four nodes operate in half-duplex mode, 
and that the destination can not communicate with the source directly and the two relay nodes 
will not interfere with each other [7]. Now we first review the two cooperative modes, namely 
SRP and AFP, as follow. 

SRP: A time block, T, is divided into two stages. 

1) Stage T. In the first AT slots, source node N s and relay node N 2 transmit signals while 
relay node N± and destination node N d are in receiving status. A G [0, 1] is time sharing 
parameter. 

2) Stage 2: In the remaining (1 — A)T slots, source node N s and relay node Ni transmit 
signals while relay node N 2 and destination node Nd are in receiving status. 

AFP: A time block is also divided into two stages. 

1) Stage T. In the first AT slots, source node N s transmit signals while both relays Ni and 
N 2 are in receiving status. In this case, the destination can not hear the signal. 
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2) Stage 2: In the next (1 — A)T slots, both relays Ni and N 2 forward the received signals in 
the first AT time slots while the destination is in receiving status and source node keeps 
silent. 

Let Xi(t) and Yj(t) denote the signal sent by node % and that received by node j, respectively. 
Then the equivalent baseband signals of the two cooperative modes are given by, respectively, 

1) In SRP: 

In the first AT slots: Y 1 (t) = g sl X s (t) + Z sl (t) , Y d (t) = g 2d X 2 (t) + Z 2d (t) 

In the remained (1 - A)T slots: Y 2 (t) = g s2 X s (t) + Z s2 {t) , Y d (t) = g^X^t) + Z ld (t) 

2) In AFP: 

In the first AT slots: Y x (t) = g sl X s (t) + Z sl (t) , Y 2 {t) = g s2 X s (t) + Z s2 {t) 

In the remained (1 - A)T slots: Y d (t) = g ld X x {t) + g 2d X 2 (t) + Z 12d (t) 
where {Z s i(t), Z s2 (t), Zi d (t), Z 2d (t), Zi 2d (t)} are independent and identically distributed (i.i.d) 
circular Gaussian random variables, is the gain factor of the link from node i to node j. In the 
discussion, all the nodes i £ {N s , Ni,N 2 , N d } are power limited. Their maximum transmitted are 
assumed to be the same, denoted as P c . Here the channel side information is used by receivers 
feedback few bits reflecting the link conditions. The bit number is relatively little compared to 
the data packets. For instance, one relay can first estimate the status of the two links associated 
to it and then feedback to the source and the destination. The source will inform the status of 
this link to another relay in next time slot by adding a overhead in its traffic massage. This 
process renews only once at the beginning of each time block according to the block channel 
fading. In this way, all the nodes could obtain the channel side information and cooperate in 
this time block, T, which consists of many time slots. 

B. Finite State Fading Channel 

To effectively analyze the fading and time-variant characteristics of channels, a finite state 
fading model is built by partitioning the instantaneous received SNR into N levels. Some partition 
methods , such as the equal-probability partition method [20], the optimum Minimum Mean- 
Square Error (MMSE) Lloyd-Max quantification method etc. have been employed previously. 
For simplicity, we shall adopt an equal-probability partition and use its mean value of SNR to 
represent the exact SNRs in each interval [20] here. Note that it will be consistent with the real 
case as the partition level is large enough. 

Let ix i (1 < i < N) denote the probability of link state i and B(i) (1 < % < N — 1) denote 
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the partition boundary of each state, which is determined by 

1 

7Ti = 7T 2 = • • • = 7ryv_i = H N = — (1) 
fl(t+l) i pB(l) i p+oo i 

f(x)dx = -;ie[l,N-2], f( x )dx = -, / /(x)dx = - (2) 

B(i) iV JO iV JB(N-l) iV 

where /(x) = Ae~ Az is the p.d.f. of the received SNR over a Rayleigh channel and A is the 
reciprocal of the average received SNR. The mean value of SNRs in each interval, SNRi,(l < 
i < N — 1), is given by 



Jl^x-f{x)dx = n rB(i + D 
Q +1) f(x)dx 4) 
The maximum average rate Cj of the link at the state level i is given by Cj = \ log(l + SNRi 



f B( ' +1) x . f( x \ d x rB(i+l) poo 

SNRi = il L+i) = N x f( x "> dx ' SNRn = N dx (3) 

L(i) f(x)dx JB(i) Jb(N-I) 



2 

In the paper, we assume that the maximum rate can be approximately achieved by some 
effective pseudo-random channel coding schemes and that the destination node can receive the 
signals from the two relays coherently in the AFP model. Thus, we use the corresponding 
capacity to approximately trace the maximum achievable transmitting rate. For simplicity, let C 
denote the maximum achievable rate in the sequel. 



C. Block Rayleigh Fading Model 

It is assumed that the channels endure block fading, which means that the received SNR in 
one time block T is a constant, but it may vary from block to block obeying the exponential 
distribution, which is corresponding to Rayleigh fading. Furthermore, the states of the four links 
are assumed to be independent and for each link its SNR varies according to a memoryless 
mode. 



D. Marshall's Queueing Theory 

For G/G/l model, Marshall's theorem on the estimation of the average queuing time is reviewed 
here [11]. 

Theorem 2.1: For all G/G/l queues with p < 1, we have 

1 ; 2A(l-p) 2v h 

where a and b denote the arrival interval and the service time, respectively. o 2 a and of denote 
the corresponding variances of them, respectively. A is the average traffic arrival rate and p is 
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(2) 

the traffic intensity of the system. Vh and v h are the first and second order moments of the idle 
period h of the system. 

v (2> 

If the inequality > j(l — p) holds, the upper bound of the waiting time is given by 

2A(W) 2A (1 " P) - 2(1 -p) (5) 

Note that the inequality (5) becomes an equality when p approaches 1. 



III. Opportunistic Scheduling Scheme Combining SRP and AFP 

To compare the performance of cooperative modes, SRP and AFP, we first analyze the capacity 
of the SRP and AFP. For convenience, several symbols are defined first. 

(1) is defined as G i5 = Pc| ^ J " 2 and C tj = \ log(l + Gy) 

(2) x,y E M + are defined as follows, respectively, x = CidC 2 d — C s \C s2 , y = C s2 C\d — C s \C 2 d. 



A. Capacity of SRP and AFP Modes 

The SRP mode is a 2-hop strategy in which relay nodes decode their recdeived information 
first before re-transmitting to the destination. It is an efficient cooperative scheme for the diamond 
relay model due to the full spatial reuse. In the AFP mode, both relay nodes just amplify the 
signals received in the first half of T and re-transmit it in the next one. At each relay node, 
the signal is multiplied with a constant and the amplified signals from the two relay nodes 
are coherently added up at the destination if the timing synchronization and carrier recovery 
are perfect. In both modes, no buffer is used by the relay nodes. The maximum achievable 
transmission rate is based on the capacities between the links associated to the relay nodes. 

Theorem 3.1:[7,Theorem4.1 and 5.2] (i) In SRP mode, (Cy, % G {s,l,2},j e {l,2,d}) 
denotes the capacity of the link from node i to node j. Transmitting rate of the link N s — N\ — Nd 
is denoted as C\ and the one of the link N s — N 2 — Nd is denoted as C 2 . The capacity of the 
SRP mode denoted as Csr, which also represents the maximum achievable rate, is given by 

Csr = max{Gi + C 2 } = max{(AiG s i + min{AiG 2d , (1 - Ai)G a2 }), 

(A 2 G 2d + min{A 2 G ld , (1 - A 2 )G sl })} (6) 

A2 

in which Ai = C ld /(C sl + C ld ) and A 2 = C s2 /{C s2 + C 2d ). 

In addition, the link-state space can be divided into four different subspace according to the 
following conditions (7) (8) (9) and (10). The explicit expression of capacity for the SRP mode 
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in each case is given as follow. 

// (x > f]y > P| yd > xC s2 ) \J(x > P| y < p| \y\C ld > xC s2 ) 
C sl (Ci d + c s2 ) 

^^sr- — ~ — V) 

// (y > p| x > P| xC s2 > yd) \J(V > pi x < P) |x|C M > yC 8 i) 

. n _ c 2d (c sl + C s2 ) 

U SR - (o) 

<~<2d + ^s2 

If (x < P y > p |x|C M < yC sl ) \J(x < P y < f] \y\C u > \x\C 2d ) 

=^ = — n r (?) 

J/ (y < P x > P \y\C ld < xC s2 ) \J(y < f] x < f] \y\C ld < \x\C 2d ) 

. ^ C , s2 (Cid + C 2(i ) 
=^ = — n r (10) 

^2d + tvs2 

(7z) In AFP mode, parameter a and /3 denote the amplified factors at relay node N\ and N 2 , 
respectively. Since the signals are received coherently, the maximum achievable rate is 

C AF = max{- • -log(l + , , o 2 , 1 -g)} 

a,/3 2 2 or + p z + 1 (7 

< max{-.-log(l + a 2 + p 2 + 1 -2)} (ID 

The first factor | is due to the equal time-sharing and the first inequality becomes equality when 
a//3 = I \g s x\ \/\\g S 2\ \ holds. 

Based on the theorem above, one can compare the maximum achievable rate of the two 
cooperative modes. However, Theorem 3.1 only consider the case where all the link capacity 
are fixed. If all the links are time varying, it is possible for us to select an effective processing 
mode adapted to the variation of the links so that we can get larger throughput by using buffers 
at relay nodes, which will shown later. Numerical results in Fig. 5 indicate that the upper bound 
of Caf in Eqn.(ll) is much smaller than Csr on average. In addition, a general form of Caf 
will be given in Appendix (B). 



B. Comparison Between SRP and AFP and An Opportunistic Scheduling Scheme 

Using Eqn.(6) and Eqn.(ll), we can divide the link-state space spanned by the four channel 
gain factors, {g s i, g S 2, 9id, 92d], into eight different subspace. In each subspace, the capacity for 
both SRP and AFP modes are completely determined. Now, let us see a special case on AFP. 
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Theorem 3.2: In AFP mode, the link-state space is divided into two different subspace by 
the following conditions, where a/j3 = \\g s i\\/\\9s2\\ is satisfied. The corresponding explicit 
expression of Caf in each case is given by 
j j. C2d ^ G\d 



Caf =---log(l + — •- — — ^ r - Ti — rT ^— ^ — — — - (12) 



G S 2(G s2 + l) " G s i(G s i + l 

1 1, , Pr 

2"2 l0g(1 + a 





9id\\ 


2 (ll^l 


ll 2 + 


\9s2 


2 ) 2 




l 2 (ll^l 


l 2 + 


11^2 1| 


2 ) + ll 


9si\ 


2 (IM 


12 i Si) 
1 ^ Pj 



G, 2 (G, 2 +1) G„(G sl + l 





^m|| 


2 (ll^l 


ll 2 + 


\9s2 


2 ) 2 




l 2 (ll^i 


l 2 + 


llfl'aall 


2 ) + ll 


9<a\\ 


2 (ll^2| 


12 , 2l) 
\ ^ P c ) 



^AF = - • - log(l + — ■ 7 ^ 77x 77^— r, . ~ ^777, ,,n . ( 13 ) 



The proof is given in Appendix (A), while a general expression of Caf is given in Appendix 
(B). 

Based on Theorem 3.1 and Theorem 3.2, the maximum average achievable rate can be 
explicitly presented in theory for all the different subspace so that the comparison between them 
becomes traceable. 

The whole link state space spanned by {g s i, g S 2, gid, g2d} is divided into eight different sub- 
space and they are given by 

subset(a) = {g s i, g s2 , gid, 92d | (7)&(12) hold.}; subset(b) = {g s i, g s2 , gid, g2d | (8)&(12) hold.} 
subset(c) = {g s i, g S 2, gid, g2d \ (9)&(12) hold.}; subset(d) = {g s i, g S 2, gid, g2d 

| (10)&(12) hold.} 

subset(e) = {g s i, g S 2, gid, g2d \ (7)&(13) hold.}] subset(f) = {g s \,g S 2,gu,g2d | (8)&(13) hold.} 

subset(g) = {g sl , g s2 , gid, g2d I (9)&(13) hold.}] subset(h) = {g sU g s2 , gu, g2d | (10)&(13) hold.} 

According to the above analysis, an effective hybrid scheme combining the SRP and AFP 
modes is presented here. 

An Opportunistic Scheduling Scheme : In the diamond relay network model, if the channel 
side information is obtained by all the nodes at the beginning of each time block, it is possible 
to find the thresholds properly for subspace partition and select an effective cooperative mode 
with larger throughput. 

Note that the switch between the two cooperative models will not be very frequent due to the 
block fading channel. It only occurs between two time blocks if necessary. 
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IV. A New Processing Strategy With Buffers Employed by Relay Nodes 

Numerical results in Fig. 5 indicated that SRP always performs much better than AFP on 
average due to its full spatial reuse. This means that only adopting a hybrid scheme simply can 
not bring an expected throughput improvement. Since the links associated to the relay is time- 
varying and block fading, some links in bad conditions will form a bottleneck, which greatly 
reduce the throughput. To overcome it, we shall put forward a new processing strategy in which 
buffers are used at both relays to efficiently use the links in good conditions, which is referred 
to opportunistic scheduling. 

In fact, it is possible to find an effective opportunity scheduling under some cases, such as 
that only one or two links are in very good conditions while the others are all in bad conditions. 
Now we consider the following two cases shown in Fig. 3. 

1) At least one link from N s to the two relays are in pretty good conditions while the two 
links from the relays to destination are in rather bad conditions. 

2) At least one link from the two relays to destination are in pretty good conditions while 
the two links from source to the relays are in rather bad ones. 

In Fig. 3, symbol G denotes the link in pretty good conditions, B denotes it in rather bad 
ones and X denotes it in arbitary conditions (i.e. pretty good, rather bad or average). In this 
two cases, the links in rather bad conditions form the bottleneck of the networks, especially for 
those schemes without buffers. Therefore, the new opportunity scheduling scheme is expressed 
as follows. 

New Strategy With Buffers Employed 

(1) Under the above tow cases, we shall use these links with pretty good conditions to transmit 
signals while the links with rather bad conditions will keep silent and some received signals will 
be stored at its corresponding buffers. 

(2) For other cases, we still adopt the same policy as that in hybrid scheme of combining 
SRP and AFP without buffers. 

From the new opportunity scheduling scheme, one can find that a new raising problem is how 
to determine each link condition being good or bad? Furthermore, the evaluation of queuing 
delay caused by using buffers is a new problem to be considered. In the sequel, we shall deal 
with them. 

A. Two Criteria for Selecting Threshold of Link Condition 

Since the signals transmitted from source to both relays in AFP mode is the same and the 
destination needs to receive it coherently, no buffers in AFP mode may loss some opportunities 
to adapt the link time- varying conditions. For SRP mode, source will transmit different signals 
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to the two relays, the sequence number of the received symbols at the destination may be 
different with its original ones due to these link time varying capacities, which is similar to the 
phenomenon in Internet. Furthermore, Csr is mainly determined by the very bad links according 
to Theorem 3.2 when the conditions of the four links have much difference. Therefore, buffers 
employed in SRP mode for the two cases listed above can mitigate the impact of bottleneck 
links. In other cases, e.g., the conditions of the links are almost the same, no matter whether it 
is good or bad, the impact of bottleneck links is relatively small for SRP mode. 

To separate the pretty good and rather bad link states, we first propose two thresholds of the 
link state levels, denoted as C upT H and C dwT H, respectively, based on the finite state fading 
channel model. That is to say, the link is considered to be in pretty good condition if its state 
level is above or equal to C up th and considered to be in rather bad condition if its state level 
is below or equal to C dwTH . 

Criterion 1: C upTH > 2 ■ C dwTH (14) 

Based on this criterion, we shall prove that only under the two cases discussed above, the 
throughput improvement can be obtained. The detail of proof is given in Appendix (C). Here 
we shall give an explanation for Criterion 1 in principle. 

Consider two consecutive time blocks T\ and T 2 , T\ = T 2 = T, and T\ is prior to T 2 . 
For the hybrid scheme without buffers, the total amount of information transmitted in this two 
consecutive time blocks, denoted as THR\, is given by 

THR l = T x • max{C 5 K(r 1 ), C AF (T X )} + T 2 • max{C SR (T 2 ), C AF (T 2 )} (15) 

and for the new developed strategy, the total amount of information transmitted in this two 
consecutive time blocks, denoted as THR 2 , is given by 

THR 2 = T ■ max{mm{C , sl (T 1 ), C ld (T 2 )}, min{C , s2 (T 1 ), C 2d (T 2 )}} (16) 

Note that "mm{C S i(Ti),Ci d (T 2 )}, {i = 1,2}" is to guarantee that the buffer is not empty 
while "max" in Eqn. (16) means utilizing the better route. In fact, in a causal system, one can 
not obtain C id (T 2 ) during the block T 1? thus Eqn. (16) only presents an ideal case. In practice, 
we shall use the better front-side link associated to the two relays when case (1) happens. This is 
because the probabilities of different link states are equal and the link states change independently 
from one block to next, which means the two back-side links associated to the relays have the 
same probability being in pretty good conditions when case (2) happens. From statistic view of 
point, choosing the better front-side link and sending massage as much as possible will obtain 
a larger gain. In addition, if the buffer becomes empty as case (2) happens, the relay will not 
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transmit any massage. This is the case just as that in a general G/G/l system. Even though, 
the performance will degrade little because the throughput of those non-buffer strategies also 
become rather small when case (2) happens. Thus, we can still use Eqn.(16) in the following 
analysis approximately. 

Let us observe the case that in T\ only the link from N s to N\ is pretty good and in T 2 
only link from N± to N d is pretty good while all the other three links are rather bad, which 
means C si (Ti) > C up th and Cid(T 2 ) > C up th hold. In this case, the throughput of the hybrid 
scheme without buffers is very small due to the existence of bottleneck links. But for the new 
developed strategy, since one buffer is employed by Ni, the total throughput in the two time 
blocks becomes much larger due to the efficient utilization of the link from N s to N\ in T\ and 
the good link from N\ to Nd in T 2 . 

Associated with that SRP mode is better than AFP on average, it can be concluded that if 
the inequality THR 2 > THRi holds, then the total amount of information transmitted in these 
two consecutive time blocks with the buffer's help will be larger than that in the original hybrid 
mode in a certain degree. For other cases except cases (1) and (2), buffers will not be used 
and their corresponding throughput parts keep the same. Therefore, one can see that the new 
developed strategy will improve the average throughput. The proof of THR 2 > THRi is given 
in Appendix (C). Its main idea is that Csr < CdwTH holds in these two cases according to 
Theorem 3.1. If the condition C up th > 2 • CdwTH m Criterion 1 is satisfied, the inequality 
THR 2 > THRi will hold an d the improvement of throughput can be guaranteed. 

TABLE I shows the levels of 16- state partition of the fading channel obtained by the cor- 
responding maximum achievable rate with the normalization P c /cr 2 = 1 for different received 
SNR. According to Criterion 1 , CdwTH is selected from the left side of * for each link state level 
and C up TH is selected from the right side of * for the corresponding state level. In addition, the 
two sides are formed symmetrically because the performance of original schemes is relatively 
good for the case that link state levels fall into the interval between * and * since the impact of 
the bottleneck only dominates when the link states differ a lot. 

B. Delay Analysis 

In the previous example, the two time blocks, T\ and T 2 , are assumed to be consecutive. In 
fact, they may not be adjacent to each other according to the i.i.d link state model. That is, the 
information in the buffer has to wait for transmission. On the other hand, due to the links being 
in Rayleigh fading, the probability of links in pretty good conditions is relatively low, resulting 
in a larger delay for information transmission. In the new developed strategy, the mean value 
of the delay for information transmission is mainly determined by the thresholds C up th and 
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CdwTH- Based on the topological symmetry of the diamond relay networks, we only analyze the 
performance of the subsystem, shown in Fig. 4. In Table II, some notations are firstly defined, 
with the i.i.d link state model, P x and P y are given by 

N -U+l.U -1 1 N — U + 1.. d 2 N 2 -(U-lf f d^ 
p — < 1 ) — y — : '—{ — Y 

N N 2 N JK N J 2N 2 K N J 

N-u + 1 ,u-\ liV-u + l w D N9 N 2 - (u-1) 2 .D^ 
P = ( 1 )( — ) = — ( — r 

y N N 2 N JK N J 2N 2 K N J 

Note that the item N ~E +1 ( \ N ~^- +1 ) is the probability that the condition of the link between N s 
and i?i is better than that between N s and simR 2 when they are both pretty good. In the new 
developed strategy, only if the link state falls into one of the two cases mentioned previously, 
the buffer works. In other cases, the buffer will not work and the networks run with the same 
procedure as the non-buffer scheme. That is, no new massage is put into the buffer though it 
is delivered to the corresponding relay node. The massage queuing in the buffer previously will 
wait for the moment at which cases (1) or (2) happen again. The details of the queuing model 
can be described as follows. 

1) Arrival process: Since the source always has information to deliver, the link state to the 
relay determines the input process of the buffer, refer to the arrival process of buffer, 
including the arrival interval and rate. When link state is at level k, the arrival interval 
between the successive traffic units is equal to e/Ck- The symbol e E (0, 1] is a parameter 
determined by the traffic types of source (i.e. bit, byte or packet), which is not a key point 
and for simplicity, let e = 1 represent a packet. 

2) Service process: The service process refers to the delivering process of the message stored 
in the buffer of the relays to destination. Once a packet stored in the buffer is transmitted 
successfully, it is served. Thus, the link state between the relay and destination determines 
the service process, including the service interval and rate. 

The arrival interval of traffic units needs to be considered for the following three cases: when 
the buffer works, the condition of the link to the relay is pretty good. When the buffer works, 
the condition of the link to the relay is rather bad; and that buffer does not work. The value of 
arrival interval and its corresponding probability for each case is 

' 1/C h prob. equal toP x /(N -U + l), iE [U, N] 
< 1/Cj, prob. equal toP y /D, j E [1,1?] 

nT, prob. equal to (1 — P) n P, n E [1, +oo), n E Z + 
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and for the service time, similar results are given by 

l/C h prob. equal toP y /(N - u+ 1), i £ [u, N] 

l/Cj, prob. equal toP x /d, j £ 

nT, prob. equal to (1 — P) n P, n £ [1, +oo), n £ Z + 
It is easy to see that the distributions of the arrival and service interval do not obey the uniform 
or Poisson distribution,etc. They are general. Therefore, one can solve it with the G/G/l queuing 
model. The mean values and variances are given in Theorem 4.1. 

In fact, the stability of queueing in buffer is very important. To solve this problem, the traffic 
intensity of buffer, denoted as p, should be less than 1. In the new strategy, the stability of buffer 
can be guaranteed by selecting the thresholds, C up th and CdwTH, following criterion 2: 

Criterion 2: (i) U = u and D > d ; (ii) U > u and D = d (17) 

The proof is given in the Appendix (D). Furthermore, one can easily infer that for the case 
"U > u and D > d" or "U > u and D > cP, the buffer is also stable. Besides, due to that 
the thresholds, U and u or D and d are relatively close to each other and the source node has 
a continuous traffic to deliver, the arrival rate of the buffer is then close to but smaller than the 
service rate. That is, the traffic intensity of the buffer, the ratio of the arrival rate to the service 
rate, approach to 1. Consequently, the upper bound of the average delay in Marshall's theory 
become more effective. 

Theorem 4.1: For the new developed strategy with buffers at both relay nodes, we have 

E(a) = E(tu)P x + E{t D )P y + ((1 - P)/P)T, E{b) = E(t u )P y + E(t d )P x + ((1 - P)/P)T(18) 
a 2 a = (E(t 2 u ) + E 2 (a)-2E(a)E(t u ))P x + (E(t 2 D ) + E 2 (a)-2E(a)E(t D ))P y 



-(1 - P)E 2 (a) + (^)2T ■ E(a) + ( (1 P) £ P) )T 2 (19) 



°l = (E(t 2 u )+E 2 (b)-2E(b)E(t u ))P y +(E(t 2 d )+E 2 (b)-2E(b)E(t d )) 2 P 



X 



+ (1 - P)E 2 (b) + ( l -j^)2T ■ E(b) + ( (1 P) )T 2 (20) 

w < — °1±A — on 

- 2(E(a) - E(b)) { ) 

where E(tu) , E(t D ) , E(t u ) and E(td) are given by, respectively, 

N D N d 

£i/C, £i/C, £i/C, £i/C, 
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The proof is given in the Appendix (E). Theorem 4.1, presented the upper bound of the average 
delay, W, which is a function of the thresholds U , D, u and d. Consequently, it provides 
an theoretical way to consider the good tradeoff between the average delay and the network 
throughput. 

C. Tradeoff Between the Throughput and the Delay 

Follow the above discussions in Subsections IV. A and IV. B, the new developed strategy can 
really improve the throughput by efficiently utilizing the pretty good links at a cost of queuing 
delay. Using the two criteria, for the fixed D and d, if a relatively lower u and higher U are 
selected, the throughput improvement is less than that with selection of a relatively higher u and 
lower U. This is because that in the former case the value of (P x P y ) 2 (THR 2 — THRi) is less 
than the one in the latter case (see Appendix (F)). On the other hand, according to Theorem 4.1, 
the former case has a shorter delay due to its higher service rate. A similar result can be observed 
for selecting D and d under the condition that U and u are fixed. Simulations in Section |V] will 
also confirm this phenomenon. 

To achieve a good tradeoff between the throughput and the delay, we shall present an algorithm 
to select the thresholds, C up th and CdwTH, which is summarized as follows. 

Enumerative Algorithm '.Due to that the nodes are able to know the partition of link states 
in advance and obtain the average received SNR at the beginning of each time block T. In other 
words, they already have the side information in Table I. 

1) Step 1: Enumerate the combinations of C up th and CdwTH according to both of the criteria. 

2) Step 2: Estimate the average delay for each case via Theorem 4.1 and find out all the 
possible combinations of C up th and CdwTH whose average delay is shorter than the 
requirement by the service traffic. Let us denote the available set as T. 

3) Step 3: For all the possible pair of C upT H and C dw TH belonging to V, to achieve larger 
throughput the nodes select U as low as possible and u as high as possible or select D as 
low as possible and d as high as possible. 

V. Numerical Results 

In this section, simulation is used to demonstrate our theoretical results. In the simulation 
part, the stream traffic model is employed and the source node is assumed to have massage 
to deliver always. The average received SNR of these two relays and the destination varies 
from OdB to lOdB. The simulation period consists of 10 5 ~ 10 6 time blocks, denoted by M, 
where each time block T is equal to 1ms. For the new developed strategy, the total amount of 
information received successfully by the destination is R suc . Then, the average throughput of 
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the new developed strategy is evaluated by R suc /M. Since there are several selections of C upT H 
and C dwT H for each average received SNR, we estimate the average throughput for each case 
under the same average received SNR. For the original schemes, similar procedure are done. 

Fig. 5 shows the maximum achievable rates of SRP, AFR the hybrid scheme and the new 
developed strategy with buffers for each average received SNR. Both the theoretical and sim- 
ulation results indicate that if all the links have the same average SNR, SRP mode is always 
better than AFP mode. Consequently, the hybrid scheme simply combining both of them is 
almost equivalent to that only adopting SRP mode. Compared with the original schemes, the 
new developed strategy with buffers really improve the average throughput, e.g. there is an 
approximate increment of 0.071 unit/s at the average received SNR of 4dB. That is about 11% 
improvement. In addition, the improvement rate will decrease as the received SNR increases 
further. It is because that when the average capacity of each link increases under a higher SNR, 
the degree of the bottleneck link influence becomes smaller. 

Fig. 6 illustrates the average delay under different received SNR, evaluated in terms of the 
number of T. It is shown that for a fixed average received SNR, the average delay increases as 
the traffic intensity becomes larger and for a fixed traffic intensity, the average delay decreases 
as the average received SNR increases. For instance, when traffic intensity is 0.98, the delay 
in the average received SNR of 2dB is 69T, while the delay in the average received SNR of 
lOdB, it is only 19T. Also it can be seen that the simulated curve becomes much closer to the 
theoretical one as p approaches to 1, which demonstrated that the upper bound derived in theory 
is effective. 

Fig. 7 shows the tradeoff between the improvement of network throughput and the average 
delay, which is consistent with our theoretical predication. Some simulation results are also listed 
in Table III where the average received SNR is equal to 6dB. It is shown that if the required 
average delay by the traffic is no more than 20 time blocks, one can only obtain a throughput 
increment of 0.031 lbit/s when the thresholds are selected as D = 3, d = 2, U =16 and u = 15. 
That is about 4% improvement compared to the original value in Table I. A larger improvement 
can be obtained if the requirement of average delay becomes loose. One can see that when nodes 
select a lower U and a higher u or select a lower D and a higher d, the throughput improvement 
become larger while the delay becomes larger simultaneously. In addition, if the average delay 
approaches to infinity in the case that D = 2, d = 2, U = 15 and u = 15, since it does not 
match the stable conditions, Criterion 2. 

VI. Conclusions 

In this paper, the diamond relay mode was studied. We compared the throughput performance 
of the two classical cooperative modes, SRP (Spatial Reused Pattern) and AFP (Amplify For- 
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warding Pattern) under wireless scenarios with assumption that the channel side information is 
known to all the nodes. We analyzed the possibility to improve the throughput by employing 
buffers at relays and proposed an new opportunity scheduling scheme. In order to improve the 
network throughput while guaranteing the stable running of the whole network, we established 
two criteria on the selection of SNR thresholds and one adjustment algorithm on the tradeoff 
between the throughput improvement and the queueing delay. Simulation results confirmed the 
effectiveness of our theoretical analysis and our new developed opportunity scheduling method. 
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Appendix (A) 

Proof of Theorem 3.2: The original expression can be transformed into the equivalent form: 



C AF = my{-.-bg(l+ 1 - 2 )} 

a 2 + (3 2 



i 2 ^ /i \\9sl \\Pc\ it \\9ld \\Pc-, a l ~ / 1 \\9s2 \\P' c\ // \\92d -P c\ 2 / n2 m l|2/|i i|2 

s.t. a < 1 + — / — , p < 1 + 2 — / 2 — h a /P = \\9si \\ \\9s2 \\ ■ 

o z a z o z a z 

That means Caf achieves its maximum when a 2 + (3 2 reaches its maximum under the three 

constraints above. With the three constraints, one can find that 

a 2 + p 2 = (1 + M| )P 2 < min{ 

. f G2d{G s i + G S 2) Gid{G s i + G s2 ) 
= GAG^TT) ' G^TTT) } <22) 



\\92d\ 


I 2 ! 


[\\9si\ 


l 2 + l 


\9s2\ 


I 2 ) 




\9s2\ 


\ 2 (\\9s2\ 


12 1 El) 

1 T pj 



\\9\d\ 


I 2 l 


[\\9si\ 


\ 2 + \ 


\9s2\ 


I 2 ) 




\9s\\ 


\ 2 (\\9si\ 


12 , <d) 
^ pj 



Thus, the expression of Caf in Eqn.(22) is obtained. □ 



Appendix (B) 



r ;1 1, , Hlffsill+IIMI) 2 ^ 

° AF = max{-.-lo g (l + a2 + (]2 + 1 ~ 2 )} 



+ 2n H^lll-PcN . \\gid \\P C a2(s 9s2 Pc \ ^ \\9ld \\Pc 

s.t. a(l + — < — , + - 2 — )< 5 — 

o z o z o z o z 
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x = a,y = (3, A = \ \g sl \\,B = \ \g s2 \\C = J — 2 — ,D 
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The original problem above is equivalent to the following one, 

(x pt, Vopt) = argmaxF, s.t.O < x < C; < y < D; A, B, C, D > (23) 

where 

F _ (Ax + By) 2 
x 2 + y 2 + 1 

To derive the maximum value of F, we first have 

OF _ 2{Ax + By) x(Ax + By) OF _ 2(Ax + By) y(Ax + By) 
dx x 2 + y 2 + 1 [ x 2 + y 2 + 1 J ' dy x 2 + y 2 + 1 [ x 2 + y 2 + 1 1 1 ' 
then we solve the following equations, 

OF A , Al+w 2 , 

_ = 0=^A + V = 5z!/=^z = - - < — (25) 

ox B y B 

d J- = ^ B + BS = A X y^ y = ^ 1 -±f)< 2 -§ (26) 

Rewriting Eqn. (23), it is easily found that 

A.l + y 2 . dF A.l + y 2 ^ OF 

0< x< - ^ — > 0, x > - ^ — < 

By ox B y ox 

This indicates that F is a monotone increasing function of x if x G (0, and a monotone 

decreasing function of x if x > ^(^-)- Similar results can be derived for 

In addition, we find that the two equations = and = can not hold simultaneously. 
Otherwise, B 2 = — A 2 (I + which resulting in contradiction. In fact, if ^ = holds, then 

\ n \x;i11 hr* miQrQntppH 



According to the expression of we only need to prove 4(1 + x 2 ) — xy > holds. Since 



dx B [ 



2 



we have 

B n , 2 , B A 1+y 2 

A {1 + x) - Xy= A + B { — 



Similar result can be obtained for that case if ^f- = holds, then ^ > will be true. 

ay ox 

In addition, according to Eqn. (25) and (26), if 2A > BC holds, then x |aF =0 > C is true and 

dx 

if 2B > AD holds, then y |sf =0 > D is true. Now we can summarize different cases: 

dy 

(1) For the case that 2 A < BC: 

If j| ( 1+ p 2 ) < C , then x opt = ^(^7^-) and y opt = D. Otherwise, we divide it into two sub- 
cases: a) If BCD > A+AD 2 holds, we have x opt = C and y opt = max{ BC+ ^^ , BC-VBW-4A* y 
b) If BCD < A + AD 2 holds, we have x opt = C and y opt = D. where "max{x, ?/}" denotes 

<2 
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the larger one between x and y with the constraint that both of them are less or equal to z. If 
anyone of them is above z, its value is defined as 0. 

(2) For the case that 2B < AD : 

If ^(^tj-) < D, then we have x opt = C and y opt = -f Otherwise, we also divide it into 

two subcases, a) If ACD > B+BC 2 , then we have x opt = max{ AD+v ^ -^^, AD-VAW-w y 
and y opt = D. b) If ACD < B + BC 2 , then we have x opt = C and y opt = D. 

(3) For other cases: F achieves the maximum value if x opt = C and y opi = D. □ 

Appendix (C) 

Without loss of generality, let us consider the proof of case (a) shown in Fig. 3. ( i.e. At least 
a link from source to the relays is pretty good and the two links from the relays to destination 
are rather bad.) It is easily to check that the condition (9) in Theorem 3.1 is satisfied for case 
(a). Thus, the maximum achievable transmitting rate can be given by 

n C ld (C sl + C 2 d) , n „s 

C S R = n n U/) 

where C sl > C upT H , C ld < C dwT H and C 2 d < C dw TH holds in case (a). 

In addition, Csr is a monotone increasing function of both C\d and C 2 d- Thus, we have 

n ^ CdwTH(C s i + CdwTH) _ n . 00 . 

<^SR < 7^ — ^ — ^dwTH Uoj 

^sl + ^dwTH 

Likewise, similar proofs can be given for another three cases in Fig. 3. Finally, based on 
Criterion 1, we conclude that Csr < CdwTH holds only in the cases (1) and (2) presented in 
Section IV. □ 

Appendix (D) 

Firstly we deduce the average arrival rate and service rate for the system shown in Fig. 4, 
denoted as A and u,, respectively. According to the new developed strategy, only cases (1) or (2) 
happens, the buffers start to work. Therefore, the average transmitting rate when the link states 
are in pretty good and rather bad conditions can be derived as follow 

EN f~i C *r^d ^i 

i=U °i Y*u _ 2-^i=u °» atD _ Z^i=l °' (~id _ 2-^i=l °i 
N _ u + l ^H~ N _ u + l ^L~ D d 

Associated with the probabilities that the buffers is in working status, P x and P y , we have 
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Next we consider the two conditions expressed in Criterion 2, respectively, 
(i) If U > u and D = d, C*f = C* = ~C~ L holds, then we have 

_ ( iV M °^ 2A^ ° H 2iV^ +Gi( 2iV^ 27V 2 " jJ 

[ N } 2iV 2 2iV 2 +l 2 )[ 2iV 2 

iy 2 -(f/-l) 2 n _ /v 2 -(f/_i) 2 iv 2 -(u-l) 2 



-)] = 2(-) 2 (^ - Q)( -^L J£-^L) < 



2iV 2 yj v iV y v w ayv 2iV 2 2iV 2 
The first inequality above holds because C L < C dwTH < \C upT H < |C]j < 1S true - 

Associated with — C H > and — ~ 2 ^~ 1) — < 0, one can see that X < fi. Thus, 

the corresponding traffic intensity is less than 1, indicating that the queueing is stable, 
(ii) If U = u and D > d, Cf = C% = Ch~ holds, then we have 

Similarly, the first inequality above holds because Cf < Cf < CdwTH < \C up th < \Ch is 
true. Check all the possible combinations in Table I according to Criterion 1 and U = u, D > d 
for each SNR case, one can see that Cf (^) 2 ) - C^(-^) 2 ) < is true. 

Likewise, we can also verify the case of the 8-level and 32-level partition and make sure the 
conclusion holds based on the two Criteria. The explanation in principle is that the impact of the 
variation of the terms "_D 2 " and "<i 2 " to the plus-minus of the inequality is much more greatly 
than that exerted by the gap between the Cj's value in [Ci, CdwTH] for each case with equal 
partition. Thus, X < fi holds, which guarantees the corresponding traffic intensity is less than 1 
and the queue is also stable. □ 

Appendix (E) 

Proof of Theorem 4.1: Here just give the proof of the mean and variance of arrival interval 
a, respectively, denoted as m(a) and <r 2 . Similar proof for the service time can follow this one. 
(1) When the buffer works and the link state is higher or equal to state level U, we have 

E fa) = £ 77 ■ Estate level = i} = £ i • l ^^" 1 ^ = P* " to (3D 

i=U 1 i=U % 
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(2) When the buffer works and the link state is lower or equal to the state level D, we have 

C j ' i—-— J j 2 iV 2 iV 

(3) When the buffer does not work in the consecutive n time blocks, we have 



\ nr , , , A 1 1 N 2 — (u — l) 2 -D „ _ 



CxD 00 1 7^ 

E(tu 3 ) = ^^ n ^ " arrival in continuous nT} = E nT ' i 1 ~ P T P = ( ) T ( 33 ) 

n=l n=l 

With the results above, we get E(a) = E(t Vl ) + E{t u% ) + E(t U3 ). 

The proof for the variance of the arrival interval is similar. For example, When the buffer 
works and the link state is higher or equal to state level U, we have 

N N 

< = Efe -E{a)) 2 -P{state level = i} = {Y,{ 7 r) 2 + {N-U + l)E 2 {a)-2E{a)- 

i=U 1 i=U 1 

E ?V) • Jj N+ 2N~ l { Jj )2 = Px ' ml) + E2{a) - 2E ^ E ^ 

i=U 

When the buffer works and the link state is lower or equal to state level D, the proof of of is 
similar and we omit it here. When the buffer does not work in the consecutive n time blocks, 
we have 



^^[(nT — E(a)) 2 ■ P{no arrival in continuous nT}] (35) 

n=l 

oo oo 

E[(™ T - E ( a )) 2 ■ C 1 - p r p ] = E^ 2t2 + p2 (°) - 2raT ■ p (°)) ■ ( x - p ) np ] 



n=l n=l 

By using the known results, for q £ [0, 1), J2^=i n2 Q n = n i^p = 
and 



m 1 

E ra2 ^ = _!)3 • -g + g m + 2mg m + m 2 g m 

n=l ' 

Y l+m o^v,^, 1 + m 0^v, 2 ^l+m , ,2,,2+m 



we have 



+g 1+m - 2mq 1+m - 2m 2 q 1+m + m 2 q 2+m )], (36) 



< = (1 " ^)^(a) + (^)2T ■ ^(«) + ( (1 ^ P) )T 2 (37) 

r 2 _ „2 i —.2 i „2 



Finally, we get o^ = o^ + o^ + o^ . □ 
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Appendix (F) 

Consider the average throughput improvement in the two time blocks Ti and T 2 mentioned 
in our discussion. It can be denoted as P{buffer, works} ■ (THR 2 — THRi). In addition, 
according to the expression of THR 2 , we need to find out the transmitting rate (i.e.the explicit 
value of min{C si (Ti) , C i( i(T 2 )} ) in all the possible cases and their corresponding probabilities. 

For the fixed D and d, if we choose U £ [C upT H + 1, -W], then u £ [C upT H, U] subjected to 
Criterion 2. The probability that the buffer works, P new , in the two time blocks T\ and T 2 is 



_ N 2 -(U-lf d N 2 — (u— I) 2 D 2 

In addition, if C si (Ti) > C id (T 2 ) holds, the average transmitting rate, denoted as A, is 

t/-i 1 N-l AT - ' 

A = Y1 Af-„4 1 C! + S (Ar-TT+lVK-„. + rt Ci (39) 



iV-u+1 ^ + l)(iV-tt + l) 



and if C si (Ti) < C ic i{T 2 ) holds, the average transmitting rate, denoted as B, is given by 



N 



Summarizing the results above, we have THR 2 = (A + B) ■ T. That is, 

1 (7-1 jv-i 1 

THR 2 = \ y^C k + ( J y(2N-2k + l)C k + C N \- -}-T (41) 

2 [ N-u + l^ *^Z^ ^ > k ^ (N — U + l)(iV — u + 1) V 7 

Let ^aj,u) represent the throughput improvement when U and u are selected. Then 

*(!/,«) = Pne W (THR 2 - TMi) = (j^Jj) 2 - ~ ~ 2 ^2 ~ ' 

(JV-t7 + l)(jV-u+l) 

Next we keep u as a invariable. For the case that U' = U + m{l<m<N — U), and we 
have 

*(tf+m,«) - Pnew{J-HR 2 - IHRi) -{ — —) — — 

1 ^tH^ 1 £wL. (2iV - 2ife + l)C fc + cw 



^T^ Ck+ (N-U + l)(N-u + l) - THR ^' T («) 



In order to prove that ^f(u,u) > ^(u+m,u), we use the following way. 
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Based on C u < C v < C u+m _i < C N _i < C N , in the expression of &(u,u)> we replace the 
former items (k = u, U — 1) with C u and the latter items C k (k = U, N — 1) with Cjj, 
then we obtain the infimum of \P m u \ as follows 

= S l N2 -( U - 1 ) 2 K N U _~ U +1 C » + 
1 



(jv - u+ i)(N ^TTT) l(jV " ^ " 17 + 2)c - + c » ] - THRl/T} ' (44) 

In the expression of ^m+m,u)j we replace the former items Ck {k = u, U + m — 1) with 
C[/ +m _i and the latter items C k (k = U + m, TV — 1) with C^-i, then we get supremum of 

*(j/+m, u ) as follows 

*5? +m ,„, - - (1/ + m - DH^t^,^, + (jy _ g _ m + 1 i)(jy _ u + i) g , 

(Jf-p-^jy tf - m + 2 ) _ 
(iV-f/-m + l)(iV-w + l) W1 17 J 

where5 = T(^fA) 2[A r 2 _ (M + i )2] . 

Now t 

we have 



■2iV 2 AT JV> 

Now by scaling $j^ u) and ^ +mM further according to C u < C v < C C7+m _ 1 < C^-i < CV, 



+ aSi + iv^TTT ' jTOI^ " Ti ^ T > (46) 

= 5-[N 2 -{U -1) 2 }(C U -THR 1 /T)=V 1 

*r* < ^■fiv 2 -(^+m-i) 2 if( [/+m ~ M + (iV ~ [/ ~ m)(iV ~^~ m + 2) 

< L iV iJm iV-u + l + (JV-E/--ro+l)(jV- u +l) 

-)C N -THR X /T\ 



(N -U -m + l)(N -u+iy 
= 8-[N 2 -(U + m-l) 2 }(C N -THR 1 /T) = V 2 (47) 

From the deductions above, the proof of ^(u,u) > ^{u+m,u) is equivalent to the proof of 
\&i > \P 2 - 

Using the result in Appendix (D), we have < THRi/T < 2Cd W TH- From Table 1, we find 
that (C u - THRi/T) > a{C N - THRi/T) always holds for all the cases of average SNR if 
a = |. Furthermore, 

VPi-^a > {a[N 2 -(U - l) 2 }- [N 2 - (U + m-l) 2 }} -5 -(C N - THRJT) (48) 
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This indicates that to prove \&i > we only need to analyze the expression a[N 2 — (U — 
l) 2 ] — [N 2 — (U + m — l) 2 ]. In fact, if the following inequality holds, 

m > ^a(U - l) 2 + (1 - a)N 2 + 1-U (49) 

the inequality a[N 2 - (U - l) 2 ] - [N 2 - (U + m - l) 2 ] > will be true. 

According to Table 1 , we obtain that the maximum value of v /«(f/-l) 2 + (l-a)iV 2 + 1 - U 
is equal to 0.342 when U = 14 for each case of SNR. Thus, inequality (49) is guaranteed for 
me [1, N — U\. Likewise, similar proof can be given in the case that U is first selected and u 
varies. It can be proved that when u decreases, the improvement of the throughput will become 
smaller. 

Thus, decreasing U and increasing u subjected to our established criteria will bring a larger 
throughput. 

In addition, one can see that the procedure of proof is related to the values in Table 1, which 
means it is influenced by the levels of link state partition. In fact, we can also verify the case 
of the 8-level and 32-level partition and make sure the conclusion also holds. 
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Fig. 1. Application Scenario Example. 
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Fig. 3. The case in which the new developed strategy with buffers works. 
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Fig. 4. The model of Ni and links marked with the thresholds. 




Fig. 5. The comparison of the throughput performance between different schemes. 
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Fig. 6. The mean value of delay for different received SNR (T) 




Fig. 7. The tradeoff between the throughput improvement and the average delay for different received SNR. 
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TABLE I 



MAXIMUM LINK-RATES (UNIT/S) FOR DIFFERENT SNR UNDER 16-STATE LEVELS 



SNR(dB) 



rankl rank2 rank3 rank4 rank5 rank6 rank7 rank8 





2 
4 
6 
8 

10 



0.025 0.07 0.115 * 0.16 0.205 0.255 0.305 0.355 

0.035 0.105 0.17 * 0.24 0.32 0.37 0.435 0.50 

P.0505 0.16 0.255 * 0.35 0.435 0.52 0.605 0.685 

0.085 0.24 0.375 * 0.495 0.605 0.71 0.81 0.91 

0.13 0.35 0.525 * 0.68 0.815 0.935 1.05 1.16 

0.20 0.495 0.715 * 0.90 1.055 1.19 1.32 1.435 



SNR(dB) 



rank9 ranklO rankl 1 rankl2rankl3rankl4rankl5rankl6 





2 
4 
6 
8 

10 



0.405 0.465 0.525 0.59 0.67 *0.76 0.88 1.015 

0.57 0.64 0.715 0.795 0.885 *0.99 1.13 1.28 

0.77 0.855 0.94 1.035 1.135 *1.255 1.40 1.565 

1.00 1.10 1.195 1.30 1.41 *1.535 1.695 1.865 

1.265 1.355 1.475 1.585 1.705 *1.835 2.005 2.175 

1.55 1.66 1.775 1.89 2.01 *2.15 2.32 2.495 
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TABLE II 



NOTATIONS DEFINED FOR THE ANALYSIS 



symbol 



N 
U 
a 
D 
d 
P, 



P 
a 
b 

C k 
tu 
to 

td 

E(a) 
E(b) 



P 

IT 



meanings 



total number of the state levels 
level number of C up th for source-relay links 
level number of C up th for relay-destination links 

level number of CdwTH for source-relay links 
level number of CdwTH for relay-destination links 
probability that at least one source-relay link is very good 
while the two from relay-destination links are very bad 
probability that at least one relay-destination link is very good 
while the two from source-relay links are very bad 
the sum of P x and P y 
arrival interval of the relay node 
service time of the relay node 
maximum rate of the link when its state is at level k 
set of arrival intervals when link state level is above or equal to U 
set of arrival intervals when link state level is below or equal to D 
set of service time when link state level is above or equal to u 
set of service time when link state level is below or equal to d 
mean value of arrival interval 
mean value of service time 
variance of arrival interval 

variance of service time 
traffic intensity of the buffer 
the mean value of delay 



TABLE III 

DELAY, IMPROVEMENT OF RATE, C up TH AND CdwTH FOR SNR = 6dB (16 CHANNEL STATES) 







SNR 




6dB 




average delay 


6.68 


15.08 


17.27 


29.38 


CO 


rate-improvement 


0.0199 


0.0269 


0.0311 


0.0465 




D 


3 


3 


3 


3 


2 


d 


1 


1 


2 


2 


2 


U 


16 


16 


16 


15 


15 


u 


14 


15 


15 


15 


15 
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